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In the course of its operation, the Volcano Ranch array collected data on the lateral distribution of showers 
produced by cosmic rays at energies above 10 17 eV. From these data very precise measurements of the steepness 
of the lateral distribution function, characterized by the r\ parameter, were made. The current availability of 
sophisticated hadronic interaction models has prompted a reinterpretation of the measurements. We use the 
interaction models QGSJET and SIBYLL in the AIRES Monte Carlo code to generate showers together with GEANT4 
to simulate the response of the detectors to ground particles. As part of an effort to estimate the primary 
mass composition of cosmic rays at this energy range, we present the results of our preliminary analysis of the 
distribution of r/. 



1. Introduction 

Over its lifetime, the Volcano Ranch array |Q 
collected data on the lateral distribution of air 
showers produced by cosmic rays of energies 
above 10 17 eV. This lateral distribution is param- 
eterized by a so-called lateral distribution func- 
tion. The steepness of this function, given by r\ , 
is sensitive to the depth of maximum (A max ) of 
the shower, and therefore to the primary compo- 
sition and to the character of the initial hadronic 
interactions. 

Our analysis of precise Volcano Ranch mea- 
surements of r\ is the first analysis of Volcano 
Ranch data using modern Monte Carlo tools. To 
simulate the development of the air shower we use 
the aires |U code (version 2.4.0), which contains 
the hadronic interaction generators QGSJET98 || 
and SIBYLL2.1 j|. The results of these simula- 
tions are convolved with a simulation of the de- 
tector response carried out using GEANT4 (|]. 
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2. Volcano Ranch 

The pioneering Volcano Ranch array was an 
array of scintillation counters in operation from 
1959-1974 at the MIT Volcano Ranch station lo- 
cated near Albuquerque, New Mexico. One of 
its many distinctions was the detection of the 
first cosmic ray event with an energy estimated 
at 10 20 eV ||. Over its lifetime, the array existed 
in three distinct configurations. The first config- 
uration consisted of twenty scintillator detectors 
of surface area 3.26 m 2 . Nineteen detectors were 
spaced 442 m apart. A twentieth detector was 
placed at various locations and used for the mea- 
surement of the density of particles. A second 
configuration had a larger spacing of 884 m. 

It is the third configuration which is most sig- 
nificant for our analysis. In this configuration, 
the number of detectors was quadrupled by split- 
ting up each of the twenty detectors into eighty 
detectors of surface area 0.815 m 2 spaced 147 m 
apart. This configuration, with many more de- 
tectors spaced closer together, allowed for very 
precise measurements of the lateral distribution 
of signals in the detectors. The steepness of the 
lateral distribution, and fluctuations in it, may 
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lead to an estimate of the primary composition. 

3. Simulation of detector response 

To simulate the detector response of the array 
to the ground particles, we utilized the general- 
purpose simulation toolkit GEANT4. Our proce- 
dure follows the prescription found in |1, where 
the detector response to electrons, gamma, and 
muons is simulated in the energy range 0.1 to 
1.0e5 MeV and for five bins per decade of en- 
ergy. This detector response is convolved with 
the results of AIRES air shower simulations to ob- 
tain scintillator yield in minimum ionizing parti- 
cles per square meter (mips/m 2 ). 

4. Lateral distribution function 

For Volcano Ranch data, a generalized ver- 
sion of the Nishimura-Kamata-Greisen (NKG) 
formula was used by Linslcy to describe the 
lateral distribution of particles at ground in 
mips/m 2 JTct] . This lateral distribution function 
is given as 
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Here r m is the Moliere radius, which is ~ 100 
m at the Volcano Ranch elevation. 

From a subset of showers detected by the array 
in the third configuration, the following form of r\ 
as a function of zenith angle and shower size was 
found to fit the data 0: 

N 

(r)(6,N))=ba + b 1 (sec0-l) + h 2 log 1 o(— § ) (3) 



with b Q = 3.88 ± 0.054, b x 
b 2 = 0.07 ±0.03. 



-0.64 ±0.07, and 



5. Preliminary comparison of lateral dis- 
tribution 

Here we compare the Volcano Ranch lateral dis- 
tribution with modern Monte Carlo simulations. 



In Figure [l] we show how well the simulations 
(aires/qgsjet98 with the detector response in- 
cluded) reproduce the average lateral distribution 
of 707 showers of estimated size N = 10 8 , mea- 
sured by the array in its first configuration ]l(J • 
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Figure 1. Average lateral distribution of simu- 
lated showers compared to average measurements 
from Volcano Ranch. 



As a further comparison, for fixed bins of zenith 
angle, we determine the number of particles at 
ground level from a fit to the lateral distribution 
function (with a = 1) and compare to the average 
functional form of r\ given by Equation 3. The 
results of this comparison for sec 8 = 1.1 — 1.2 
can be seen for proton and iron showers using 
aires/qgsjet98 and aires/sibyll2.1 in Fig- 
ures |2] and ||. One can see that the average form 
of 77 over a realistic range of mass composition, 
from proton to iron, is reproduced by the simula- 
tions. At this stage, any attempt to draw a firm 
conclusion about composition or the suitability of 
interaction models is premature; we have limited 
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Monte Carlo statistics and are only comparing to 
an average functional form of 77. What we are 
content to take from this stage of the analysis is 
encouragement from the fact that the simulation 
can reproduce the data. A more detailed analysis 
will follow in a later work. 
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Figure 2. Comparison of 77 as a function of shower 
size for sec 9 = 1.1 - 1.2 and aires/qgsjet98 



6. Towards an estimate of primary compo- 
sition 

A more robust estimator of primary composi- 
tion comes from an analysis of the measured fluc- 
tuations oil], which were found to be significantly 
greater than that expected from measurement er- 
ror alone 0. From our knowledge of the char- 
acteristics of showers produced by nucleons and 
nuclei, one would expect that measured fluctua- 
tions in the distribution of r\ would be smaller for 
a heavier composition. This can be understood if 
we consider an air shower produced by a nucleus 
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Figure 3. Comparison of 77 as a function of shower 
size for sec6> = 1.1 - 1.2 and AIRES/sibyll2.1 



with A nucleons as a superposition of A show- 
ers each with energy E / A. Thus in iron-initiated 
showers the average fluctuations in shower devel- 
opment are reduced. 

At Volcano Ranch, measurements of i] were 
made on a shower-to-shower basis for fixed bins 
of zenith angle g. Our procedure is to simu- 
late showers in these bins of zenith angle, find r\ 
from a fit to the lateral distribution function, and 
compare to data. The results of this preliminary 
comparison can be seen in Figures |] and |^. 

7. Conclusions and future plans 

Using modern Monte Carlo tools previously un- 
available during its operation we have been able 
to reproduce lateral distribution measurements 
taken at the Volcano Ranch array. With an 
eye towards estimating the primary composition 
from these measurements, our ability to repro- 
duce them gives us some confidence that our anal- 
ysis procedure is correct and indicates that cur- 
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Figure 4. Fitted r\ parameter for proton and iron 
compared to Volcano Ranch data with a = 1, for 
sec# = 1.0 - 1.1, using QGSJET98 



Figure 5. Fitted rj parameter for proton and iron 
compared to Volcano Ranch data with a = 1, for 
sec# = 1.4 - 1.5, using QGSJET98 



rent hadronic interaction models can describe the 
data fairly well. A preliminary analysis of the dis- 
tribution of the steepness of the lateral distribu- 
tion function, rj, indicates that indeed something 
can be said about primary composition from an 
analysis of the data, but that further analysis is 
needed to make any stronger statement. How- 
ever, it is encouraging that the data from Volcano 
Ranch may be used to say something useful about 
mass composition thirty years after its closure. 
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